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ABSTRACT 

A new technique for observer design which requires the 
solution of a minimax problem is presented. An algorithm 
for solving the minimax problem is proposed and applied to 
the design of full-order and reduced-order observers. 
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I. INTRODUCTION 



In feedback control systems, the input vector to a plant 
is often a function of the current plant state vector. This 
means that one requires knowledge of the plant state vector 
at the current time to accomplish system control. In most 
systems, however, the entire state vector is not available 
for direct measurement. Hence, it is desirable to recon- 
struct, or estimate, the unknown plant state vector. The 
device which performs this reconstruction of the plant state 
vector is called an "observer." In this thesis, after 
presenting necessary background material on the theory of 
observers, which is primarily the work of Luenberger [4], a 
new approach to observer design is systematically developed. 
This new method uses a quadratic performance measure which 
is to be minimized with respect to the observer gain matrix 
G. 

For simplicity, only linear time-invariant continuous- 
time systems are considered throughout this thesis. 

The algorithm proposed is suitable for digital computa- 
tion using appropriate subroutines and is applicable to 
higher-order systems. In order to illustrate the new 
technique, numerical examples are presented. 
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II. OBSERVER THEORY 



A. THE NEED FOR OBSERVERS 

To design a feedback control system it is often desired 
to know the entire plant state vector. But, in most cases, 
the entire state vector cannot be measured, hence a suitable 
estimate to the plant state vector must be determined. This 
creates the need for a system called an ’'observer" which 
produces an estimate of the plant state vector. 

Assuming that it is possible to build a system which can 
estimate some constant linear transformation of the system 
state vector, then, provided that the transformation is 
invertible, it is possible to reconstruct the state vector 
itself. This is the fundamental idea upon which observer 
theory is based. 

B. FULL-ORDER OBSERVERS 

First consider the construction of a full-order observer 
for a free system- described by 



x(t) = A x(t) 


(1) 


£(t) = C x(t) 


(2) 



where x(t) is the nxl state vector, y_(t) is the mxl output 
vector, A is an nxn constant matrix and C is an mxn constant 
matrix, and n > m. Thus, an m-dimensional output vector y_(t) 
is available from measurements, but we want to know the 
entire state vector x(t) of the plant. Since x.(t) cannot be 
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measured directly, one must obtain an estimate x(t), that is 
a good approximation to the original plant state x(t) for 
all time. In what follows it is assumed that the plant 
described by (1) and (2) is observable. 

One way to estimate the state vector x(t) might be to 
construct a system of the same dynamic order and structure 
as the plant, given by 

z_(t) = A z/t) (3) 

and drive it in parallel with the plant. If this is done 
the system state, z_(t), will be a linear transformation of 
the plant state x(t)* If the initial states z_(0) and x(0) 
are equal, then z^t) and x(t) remain equal for all t > 0. 
Usually, however, ^(0) is not equal to x(0) and the differ- 
ence between the two states z_(t) and x(t), can be determined 
by forming 

z(t) - x(t) = A [ z_( t ) - x(t)] (4) 

The solution of (4) is 

z_(t) - x(t) = e— t [z_( 0 ) - x_( 0 ) ] (5) 

Thus, the error z.(t) - x(t) depends on the matrix A. If the 
matrix A is a stability matrix, i.e., all the eigenvalues of 
the matrix A have negative real parts, the error approaches 
zero as time goes to infinity. 

A better alternative is to construct a linear system 
which contains a model of the plant state equation but is 
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also driven by an output error signal, y_(t) - C z^t), that 
is, 

z(t) = A z(t) + G[y_(t) - C z(t)] (6) 

or 

z(t) = (A - G C) z(t) + G y_(t) (7) 

This system has a very desirable feature: the observer gain 

matrix G can be selected to determine the rate at which the 
observer state z(t) approaches the plant state x(t). If the 
state z_(t) becomes equal to the plant state x(t) at some time 
t = t-p then the correction term, y_(t) - C z(t), is zero for 
t = t^ and the state z(t) will remain aligned with the plant 
state x(0 for t _> t-^. 

The system described by equation (6) or (7) is referred 
to as a "full-order observer" for the plant described by (1) 
and (2). The gain matrix G can be selected by the designer. 
How to select the. gain matrix G is the problem of observer 
design. 

The above development can be easily extended to forced 
systems by including the input u(t) in the observer equation 
(6) or (7)- Consider the system described by 



x(t) = A x(t) + B u(t) 


(8) 


y (t ) = C x(t) 


(9) 
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where u(t) is the rxl input vector and B is an nxr constant 
distribution matrix. The observer for this plant is similar 
to that for the free system, but the input vector u(t) must 
be included. The appropriate full-order observer system is 

z{t) = A :s(t) +Bu(t) + G[y(t) - C z_( t ) ] (10) 

or 

^(t) = (A - G C) z(t) + B u(t) + G y(t) (11) 

The configuration of this full-order observer is shown in 
Pig. 1. 




FIGURE 1. PULL-ORDER OBSERVER 



C. REDUCED-ORDER OBSERVERS 



The full-order observer possesses a certain degree of 
redundancy. The reason for this redundancy is that the full- 
order observer constructs an estimate of the entire plant 
state, but the output of the plant, which is made up of 
linear combinations of the states, is available by direct 
measurement. Hence, it is actually unnecessary to construct 
a system like the full-order observer to estimate the entire 
plant state. 

The redundancy can be eliminated by reducing the dynamic 
order of the observer system to (n - m) , the difference 
between the dimensions of the state vector x(t) and the 
output vector ^(t) . The observer which is reduced to its 
minimal dynamic order (n - m) is called a "reduced-order 
observer." The structure of a reduced-order observer is 
illustrated in Fig. 2. 




FIGURE 2. STRUCTURE OF REDUCED-ORDER OBSERVER 

The starting point for consideration of reduced-order 
observers is again the plant characterized by 



_x ( t ) = A x(t) + B u(t) 


(12) 


£(t) = C x(t) 


(13) 
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The output distribution matrix C has dimension m. By intro- 
ducing a change of coordinates, it can be assumed that the 

I 

matrix C takes the form C = [ I I 0 ] . The new state 

i 

vector 3j[(t), is selected so that its first m components are 
equal to the plant output vector y_(t) and the remaining 
(n - m) components, denoted by w(t), can be selected arbi- 
trarily but in such a way that the states y(t) and w(t) 
are linearly independent. Thus, the new state vector x(t) 
can be written in partitioned form as 



x(t) 



y(t) 
w( t ) 



and £(t) is related to the original plant state vector x(t) 
by a nonsingular matrix P, i.e.. 



x(t ) 



P x(t) 




x(t) 



The matrix C_ here is the one given in (13) and D is an 
arbitrary matrix selected in such a way that P is non- 
singular. The new state equations have the form 



1 

1 

•p 1 
♦ >i! i 

i 

1 


= A* 


y_(t) 


w(t) 




w ( t ) 



(14) 



The matrices A* and B* in (14) are not the same as A and B 
in (12). These two sets of matrices are related to one 
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another through P. Equation (14) can be rewritten in 
partitioned matrix form as 



y.(t) 

• 




A 1 ' 

-11 i 

J 


A « 

-12 




z(t) 




Sf 


w(t ) 




< — l 
<1 

1 


A 1 

—22 




w(t ) 




Sa' 



u(t) (15) 



or 



y.(t) = A 11 ' y ( t ) + A 12 ’ w(t) + B 1 ’ u(t) (16) 

w(t) = A 21 ‘ £(t) + A 22 ’ w(t) + B 2 - u(t) (17) 

The vector y_(t) is available for measurement, and if it is 
differentiated, y(t) can also be determined. Since u(t) is 
also available, ( 16 ) provides A^ 2 ' w(t) v/hich plays the role 
of a set of measurement equations for the system (17) which 
has w(t) as its state vector and A^' ^_(t) + B 2 ' u(t) as its 
input vector. The next step is to construct a full-order 
observer of order (n - m) for the plant of equation (17) 
using 

A 12 ' w(t) = y(t) - A 11 ' y ( t ) - B 1 ' u(t) 

as the measurement. Thus, the observer for the system 
described by (16) and (17) can be expressed as 
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1 



X 

w(t) 




- G A 12 ') w(t) + A 21 'y(t) + B 2 ’ u(t) 
+ G (t) - A 11 1 ^(t)] - G B x 1 u(t) 



or 

X 



w(t) 



(A 22 « - G A 12 ’ ) w + (A 21 « - G A 11 t ) y_(t) 
+ (B 2 « - G B 1 ' ) u(t) + G £(t) 



(18) 



(19) 



A 

where w represents the estimate of w. The gain matrix G 
in (18) or (19) can be selected so that the observer system 
coefficient matrix A 22 ' - G A^ 2 ' has arbitrary eigenvalues 
provided that (A^ 2 ', A 22 ') as °t>servable. It is known [4] 
that if (C, A) is observable, then (A^ 0 ’, A 22 ' ) as a - i - so 
observable. The configuration of the reduced-order observer 
is shown in Fig. 3- 




FIGURE 3. REDUCED-ORDER OBSERVER WITH DERIVATIVE 
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The differentiation of ^_(t) in (19) may be troublesome 
to perform physically; however, it can be avoided by modi- 
fying the observer state equation (19) • The resulting 
modified observer state equation can be written as 

z(t) = ( A_22 ' “ Q A 12 ' ) z(t) + ( A_22 ’ - 0 A.^' ) G J£.(t) 

+ (A 21 ' - GA 11 ') £(t) + (B 2 ' - G B x ’ ) u(t) 

( 20 ) 

with 

A 

w = z_(t) + G y_(t) (21) 

and the corresponding observer diagram is shown in Fig. 4 
which is equivalent to the observer in Fig. 3 at the point 

A 

W . 




FIGURE 4. MODIFIED REDUCED-ORDER OBSERVER 
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D. THE CONVENTIONAL METHOD OF OBSERVER DESIGN 

It has previously been shown that the full-order observer 
for the linear time-invariant plant 



x(t) = A x(t) + B u(t) 
jr(t) - C x(t) 



( 22 ) 



is characterized by 

_z(t) = (A - G C) z_(t) + G y_(t) + B u(t) (23) 

If (C, A) is observable, the observer's eigenvalue can be 
located arbitrarily by selecting the gain matrix G_. The 
problem is how to select the eigenvalue locations for the 
observer . 

The observer error e_(t) is defined as 

e(t) = z_(t) - x(t) (24) 

Differentiating equation (24), substituting the expressions 
• • • 
for z/t) and x(t) , and solving for e_(t) gives 

e_(t) = (A - G C) e_(t) (25) 

Letting 

F = A - G C 

then the error equation for the observer can be written as 
e(t) = F e(t) ( 26 ) 
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The solution to this equation is 



e(t) = e- e(0) (27) 

Thus, the dynamic response of the error between the observer 
and the plant state is determined by the matrix F = A - G C. 
It is desired that the error e(t) eventually approach zero 
for all e(0). Thus, if G is selected so that F = A - G C 
is a stability matrix then the observer state z_(t) approaches 
the plant state x(t) at a rate controlled by G. One approach 
to observer design is to select the desired eigenvalues and 
then determine G to obtain these eigenvalues [M], 
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III. A NEW METHOD FOR OBSERVER DESIGN 



A. PROBLEM DESCRIPTION 

An alternative procedure for observer design can be 
based on a performance measure for the error system (26) 
consisting of the integral of a quadratic form in the 
observer error, that is. 



where Q is a real symmetric positive definite constant matrix 
and e_(x) is a function of the gain matrix G, which is to be 
selected to minimize this performance measure. It is shown 
in the Appendix that equation (28) can also be expressed as 



where e is the initial error of the observer, e^ is the 
transpose of e^, and V^(G) is the value at t = 0 of the real 
symmetric positive definite solution V(t) of the matrix 
differential equation 



J = / e (t) Q e(t) dx 



( 28 ) 




(29) 



T 



V(t) = -F T V(t) - V(x) F - Q 



(30) 



subject to the boundary condition 



V(t f ) = 0 



(3D 
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Alternatively, defining t = t^ - t, this equation can be 
expressed as 



V(t) = F T V(t) + V(t) F + Q 



(32) 



with the initial condition 
V(0) =0 



(33) 



From equation (29) it is seen that J depends on the 
initial error of the observer, e . and on the matrix V (G) 
which is implicitly dependent on G through equation (30). 
Finding a solution to equation (30) requires solving 
n(n + l)/2 differential equations (since V is a symmetric 
matrix) . 



B. CONSIDERATION OF THE INITIAL CONDITION 

It is desired to select the observer gain matrix G to 
minimize the performance measure 



J 




e T (x) Q e(r) dx 




e 

— o 



To calculate the performance measure, a priori knowledge of 
the initial state of the error system is needed. The initial 
error state is not known, but a reasonable approach is to 
take the worst initial state in some bounded region. It will 
be assumed that the initial observer state is 0_ and that the 
possible initial error states lie within a hyperspherical 
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region with center at the origin and having a specified 
radius p. With this assumption the performance measure can 
be written as 



J 

P 



max 



( x) Q e ( x) dx 



— o 



o 



I So 



111 p 



max 







V (G) e } 
— o — -o 




P 



(3*0 



where || || denotes the Euclidean norm. It is easily shown 

that the maximum occurs where IIIqII = P > moreover, without 
loss of generality it can be assumed that p = 1 . That is, 
the maximum can be considered to occur on the unit hyper- 
sphere. This is because, for any e^ such that IIIqII = 1 and 
any radius p > 0 , 

J = p e T V (G) p e = p 2 e T Y (G) e 

p K — o — o — K — o K — o — o — — o 

= p 2 J 1 (35) 

i.e., the maximum performance measure on the hypersphere of 
radius p is p times the maximum performance measure on the 
unit hypersphere. Hence, it is desired to minimize 
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J 



1 



max 



f e T (x) Q e(x) dT 




1 



max 



e 

— o 




( 36 ) 




1 



with respect to G. 

C. EXTENSION TO REDUCED-ORDER OBSERVERS 

So far, a new method to construct a full-order observers 
has been discussed. The method can easily be applied to 
design reduced-order observers also. Consider again the 
linear time-invariant system 



The reduced-order observer equation for this system is 



x(t) = A x(t) + B u(t) 



( 37 ) 



£(t) = C x(t) 



( 38 ) 



z(t) = (A 22 ' - G A 12 ') z(t) + (A 22 ‘ - G A 12 *) G y_(t) 



+ (A 21 « - G A 11 ' ) ^(t) + (B 2 * - G B x ’ ) u(t) 



( 39 ) 
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with 



w(t) = z,(t) + G ^(t) (40) 

where w(t) is the estimate of w(t) . 

The characteristics of the observer system (39) depend 
on the coefficient matrix A 22 ' ~ G. A^ 2 ' which depends on the 
gain matrix (3. Again to determine the gain matrix G, 

A 

consider the state error, w(t) - w(t). The observer error 
e^(t) is 

e_( t ) k w(t) - w(t) (4l) 

Differentiating equation (4l) gives 

• x 

£(t) = W(t) - W(t) 

X 

Substituting w(t) and w(t) and solving for e(t) yields 

e(t) = (—22 ~ — —12*^ £- (t) (ii2) 

Letting F = A_ 22 ' “ G A^ 2 ' , then (42) can be written as 

e(t) = F e(t) (43) 

The solution to this equation is 

e(t) = c— ^ e_(0) (44) 

Notice that equation (44) has exactly the same form as 
equation (27) except for the definition of the matrix F. 
Therefore, the results obtained for full-order observers 
apply to reduced-order observers as well. 
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D. SIMPLIFICATION FOR THE INFINITE-TIME CASE 



It will now be assumed that the upper limit on the 
performance measure integral is infinity. It will be shown 
that if t f -*°° , the problem of calculating V^G) can be 
simplified to the problem of solving a set of n(n + l)/2 
linear algebraic equations. Let us state an important 
lemma and two theorems related to this simplification. 

The first step is to write an explicit solution for V(t). 

Lemma 1: The solution of the matrix differential 

equation 

V(t) = F T V(t) + V(t) F + Q ( 45) 

is 

V(t) = <j)_( t ) V(0) <ff^(t) + / £(t - t) Q ^ (t - x) dx (46) 

0 

where c£_(t) is the transition matrix for the system 

z(t) = F T .z(t) (47) 

Proof : Suppose that (46) is the solution of (45), and is 

rewritten as 

V(t) = £(t) V(0) £ T (t) + £(t) [ f £(-x) Q £ T (-x) dx] £ T ( t ) 

0 

i £(t) V(0) i T (t) + £(t) L(t) i T (t) (48) 

where L(t) is defined as 

L(t) = / t £(-x) Q i T (-x) dx (49) 

0 
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Differentiating (^8) gives 

V(t) = j>(t) V(0) ^ T (t) + ±(t) V( 0 ) i(t) + i(t) L(t) £ T (t) 

+ £(t) L ( t ) <Jt) + £(t) Ci(-t) Q i. T (-t)] £ T (t) 

(50) 

The last term of equation (50) becomes Q, because 

$(t) [±(-t) Q l T (-t)J l T (t) = 1(0) Q £ T (0) = I Q I = Q 



Also 



i(t) = F T ±(t) 



(51) 



therefore (50) can be written as 

V(t) = F T 4i(t) V(0) i T (t) + £(t) V(0) i T (t) F 

+ F T £(t) L(t) i T (t) + ±(t) L ( t ) £ T (t) F + Q 

(52) 



But from equation (^8) 

£(t) V(0) i T (t) = V(t) - £(t) L ( t ) i T (t) (53) 

Substituting (53) into (52) gives 

V(t) = F T [V(t) - £(t) L ( t ) i T (t)] + [V(t) - £(t) L ( t ) i T (t)] F 

+ F T £(t) L(t) £ T (t) + £(t) L ( t ) ^ T (t) F + Q (5*0 

Simplifying equation (5*0 yields the final desired result 
V(t) = F T V(t) + V F + Q 
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The following well-known theorem gives a necessary and 
sufficient condition for the error system (26) to be asymp- 
totically stable, that is, 

lim e_(t) = (3 for all ^(0). 

t-H» 

Theorem 1: The linear time-invariant error system 

e(t) = F e(t) (55) 

is asymptotically stable if and only if, given any positive 
definite symmetric matrix Q, there exists a positive 
definite symmetric matrix K such that 

F T K+KF=-Q (56) 

Theorem 1 has been proved in [2]. This theorem can be 
used to determine whether or not the error system (26) is 
asymptotically stable for any choice of the gain matrix G. 

The following theorem applies specifically to the case 
where t ^ ®°. 

Theorem 2: Assume that the linear time-invariant system 

e(t) = F e(t) (57) 

is asymptotically stable, and has a performance measure 
defined by 

fc f 

J = / e T (x) Q e(x) dx (58) 

0 
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which is equivalent to 



J = e T V (G) e 
— o — o — — o 



(59) 



where V(t) satisfies 

V(t) = F T V(t) + V(t) F + Q 



with 

V(0) = 0 



then 

V(t) = F T V(t) + V(t) F + Q (60) 

has a constant steady-state solution i.e., 

lim V(t) = V (a constant) 
t-H» “ 

Proof : By Lemma 1 with V(0) = 0, the solution to ( 6 0 ) is 

V(t) = / £(t - x) Q £ T (t - x) dx ( 6 1 ) 

0 

T 

Substituting <£(t) = e— ^ and <^*(t) = e— ^ into ( 6 1 ) and 
using properties of the matrix exponential, gives 



T t T T T 

F 1 t , -F x n / F 1 1 -Fx . T . 
V(t)=e— / e Q (e— e — ) dx 



F T t } -F T x D -F x . F T 
= e— / e — Q e — dx e— 



( 62 ) 
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But the system (57) is asymptotically stable, so according 
to Theorem 1 there exists a positive definite symmetric 



matrix V such that 
— o 



Q = -F V-VF 
— — — o — o — 



(63) 



Substituting ( 63 ) into (62) gives 



V(t) = - e- / e - T [F T + Vq F] e - T dx e- fc (6k) 
-F^t T 

Since e — and F commute, this can be written as 

V(t) = -e— / - 4- [e-^ V e~— T ] 

0 dT -O 



T T 

Ft r -Ft .. -Ft 
= e— Le — V e — 
— o 



n Ft 

J £- 



0 



mm m 

F x t -F t t t -Ft Ft F t T „ T Ft 

= e— e V e — e— - e— I V I e— 

o — — o — 



T *tt T F T t ,, Ft 
= I V I - e— V e— 
— — o — — o 



(65) 



The system (57) is stable, so 



, , Ft 
lim e— 

t^°° 



= 0 and 



T 

F f 

lim e- = 0 
t->-°° 



( 66 ) 



Hence, substituting (66) into ( 65 ) 



lim V(t) = V (a constant) 
— — o 

t+°° 



which proves the theorem. 
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According to the result of Theorem 2, V(t) becomes zero 

as t . Therefore. V is the solution of 

— o 

P T V+VF=-Q ( 67 ) 

which is obtained by setting V(t) = 0_ in equation (32). 
Equation ( 67 ) gives a set of n(n+l)/2 linear algebraic 
equations that can be solved for V after any positive definite 
Q has been selected for the performance measure (28). Hence, 
the problem of calculating becomes simple. The positive 
definiteness of the solution V can be checked by using 
Sylvester's criterion [ 7 ] and if and only if V is positive 
definite the system is asymptotically stable. 
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IV. SOLUTION OP THE OBSERVER DESIGN PROBLEM 



As mentioned before, the error equation needed to 
determine the gain matrix £ has the same form for both the 
full-order and reduced-order observers. So one design 
method can be applied for both cases. 



A. DEVELOPMENT OF THE METHOD 

It is desired that the gain matrix G be found which 
minimizes 






max 






I So" = 1 



m 

/ e_ (x) Qe(x) dx 
0 



max 



^o 



[ £o VS) So 3 



(68) 



HSo 



= 1 



Letting min then 



V 



mm 

G 



max 



e 

— o 



^o V G) ^o ] 



( 69 ) 



II ®r 



= l 
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This is a minimax problem, and equation ( 69 ) is a real 
quadratic form in . There is a useful theorem about the 
extremal properties of the eigenvalues of a real quadratic 
form [ 1 ], 

Theorem 3: The global maximum of a real quadratic form 

on the unit hypersphere is equal to the largest eigenvalue 
of the quadratic form, and moreover the corresponding eigen- 
vector is the vector drawn from the origin to the point on 
the hypersphere where the quadratic form achieves its 
maximum. 

Therefore 



max 

e 

— o 

Kll 



e 

— o 



Vo(G) 



^o 



] = 



X 1 (G) 



1 



(70) 



where A^(G_) is the largest eigenvalue of V^G) . Hence, now 
the problem is to minimize X^(G) with respect to the elements 
of G; that is, to find the observer gain matrix G such that 



= min X^(G) 



(71) 



To ensure that the observer is stable and does not have an 
unlimited bandwidth, the gain matrix G will be selected from 
the set JJ defined as 




1 . All eigenvalues of F have negative real parts 
G: 

2. Sum of eigenvalues of F = trace F _> -D 
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where D is a specified positive constant. For a given value 
of G_ the largest eigenvalue is easily found by using a 
digital computer. The program required is simple because 
very efficient methods are available for finding the largest 
eigenvalue of a real symmetric matrix, and many computer 
facilities have several standard subroutines that can be 
used for this purpose. 

B. ALGORITHM FOR COMPUTING THE GAIN MATRIX 

The steps in the computational procedure for finding the 
optimal gain matrix are given below. 

1. Guess G subject to the following requirements: 

a. F is stable. This can be done by assigning arbitrary 
eigenvalues with negative real parts to F and using Luenberger's 
method [4], for example. 

b . Trace F _> -D 

2. Use a minimization subroutine to minimize 



J 



1 



max 





Actually, a penalty function is added to ensure that the 



trace F > -D constraint is satisfied, e.g.. 



0 



if trace F > -D 




a x [trace F + D] 2 , if trace F <_ -D 



A 




where a is a positive weighting factor. 
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In computing values of J for the minimization subroutine 

Cl 

the following steps are carried out. 

T 

a. Solve F V + V F = - Q for V = V 

• — o 

b. Check to see that is positive definite. If it 
is continue, otherwise set J & = 10 , that is, a large 
positive number and return to the minimization subroutine. 

c. If is positive definite it is the steady-state 
solution of 

V(t) = P T V(t) + V(t) F + Q 

Next find the largest eigenvalue of using a 
suitable subroutine. 

d. Calculate the trace of F and compute the penalty 

function J . 

P 

e. Set J = J, + J and return to the minimization 

alp 

subroutine . 

In applying this algorithm to the examples given in this 
thesis the minimizations were accomplished using a subroutine 
"DIRECT" (NPS Computer Facility) that performs the pattern 
search method of Hooke and Jeeves, and the largest eigen- 
values were found using a subroutine "GIVHO" (NPS Computer 
Facility) that utilizes the Givens-Householder method. 
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V. APPLICATIONS OF THE ALGORITHM 



A. A FULL-ORDER OBSERVER 

The purpose of this example is to illustrate the appli- 
cation of the described algorithm to the design of a full- 
order observer for a second-order system. 



Y ( s ) 

Consider the plant 



S + 3S + 2 



with the state- 



variable representation 

x(t) = A x(t) + B u(t) 



— — 










0 1 




x 1 ( t ) 


+ 


0 


-2 -3 




x 2 (t) 




1 


__ __ 









- 



u(t) 



(72) 



y(t) = C x(t) 

= [ 1 - 0 ] 



^(t) 

x 2 (t) 



(73) 



A full-order observer for this system is determined by 
specifying the observer gain matrix 



G = 



fe l 

g 2 
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The observer system coefficient matrix is 



F = A - G C = 



“Si 1 

- 2 -g 2 -3 



which has the characteristic equation 



A + (3 + gqH + Sg-L + g 2 + 2 = 0 



( 74 ) 



(75) 



Selecting the observer eigenvalues arbitrarily at — 3 ± j 3 
gives the characteristic equation 



( A + 3 + j 3 ) (A + 3-33) = A 2 + 6 A + 18 = 0 



(76) 



Equating the coefficients of like powers of A in (75) and 



( 76 ) yields 



G 



3 

7 



This is the starting guessed value of the gain matrix G. If 
one arbitrarily chooses Q = I_, and solves ( 67 ) for V the 
result is 



15 + 3g 1 + 5g 9 + g 2 2 

12 + 22g x + 6g 2 + 2g x g 2 + 6g x 2 

3 ~ ~ gqS2 

12 + 22g x + 6g 2 + 2g x g 2 + 6g x 2 

3 + 3g x + g 2 + Sq 2 

12 + 22g 1 + 6g 2 + 2g x g 2 + 6g x 2 



( 77 ) 
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Due to the symmetry of V only three equations are needed. 
If the sum of the observer eigenvalues is to be limited to 
-10.0, for example, then the proposed algorithm gives the 
solution for the observer gain matrix as 



G 



7.0 

- 0.64063 



which gives the observer eigenvalues as 



( 78 ) 



X 



- 6.625 

- 3.375 



( 79 ) 



For other selections of the weighting matrix Q and different 
limiting values for the trace of F the values obtained for 
the observer gain matrix G and the corresponding eigenvalues 
are tabulated in Table 1. Looking at the table it can be 
seen that the two eigenvalues become farther apart from each 
other when the limiting value of the sum of the eigenvalues 
is decreased (assuming a fixed Q) . 

To compare this full-order observer, which will be 
referred to as the optimal observer, with other observers 
having different observer eigenvalues, consider two observers 



having eigenvalues 


selected as 






-5-0 




-8.0 


X = 




> 1 = 






-5.0 




-2.0 




, __ 




— — 



respectively. The corresponding observer gain matrices are 
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TABLE 1. OBSERVER GAINS AND THE CORRESPONDING EIGENVALUES 



G 



(81) 









— — 




7.0 




7.0 


- 




3 G ~ 






2.0 




- 7.0 




— - 




. __ 



respectively . 

For simulation purposes, the input u(t) was selected to 
be a unit step function. Assuming that the initial error 
states lie on the unit circle, the actual and the estimated 
states computed by using the optimal gains of equation ( 78 ) 
are plotted as functions of time on the graphs in Fig. 5 for 
the best initial error condition and in Fig. 6 for the worst 
initial error condition. The estimate errors for the 
observers with three different observer eigenvalues are 
plotted on one graph in Figs. 7 and 8 for the best and worst 
initial error conditions, respectively. The values of the 
performance measure for the observer error systems with 
different observer eigenvalues (different observer gains) 
are tabulated in Table 2 for the worst initial conditions. 



Observer Gains . 
(Eigenvalues ) 


V/orst Initial 
Conditions 


Values of Perform- 
ance Measure 


S 1 = 7,0 ( X 1 = " 6 ' 625 \ 
gg = -0.6i|063\A 2 = -3-375/ 


e(0) = 


0.16 

-0.987 




0.16^ 


g l = 7 '° / X 1 = " 5 '° \ 

g 2 = 2.0 ^ X 2 = -5.0 / 


e(0) = 


" 0.525' 

-O.85 




0.181 


S 1 = 7 '° I X 1 = “ 8 ’° \ 

&2 = "7-0 \^ 2 ^ -2.0 J 


e(0) = 


“ 0.63 ' 
0.777 




0.303 



TABLE 2. PERFORMANCE VALUES FOR VARIOUS OBSERVER 
EIGENVALUES 



36 





o 

M 

00 



00 

w 

EH 

< 



Eh 

00 

w 






M 

£ 

w 

Eh 



i — 1 



X 



Q 

< 



00 

W 

Eh 

< 

E-h ^ 
00 O 



^ EH 

< M 

Q 

Eh ^ 



O O 

< O 



LT\ 



ctf 



a 

o 

H 

P-c 



37 



Czti' Tfiate 




38 



FIGURE 6. ACTUAL STATES AND ESTIMATES USING THE OPTIMAL GAINS FOR THE WORST INITIAL 
CONDITION 
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FIGURE 7. ESTIMATE ERRORS FOR DIFFERENT OBSERVERS FOR THE BEST INITIAL CONDITIONS 
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(a) ESTIMATE ERROR FOR x(l) (b) ESTIMATE ERROR FOR x(2) 

FIGURE 8. ESTIMATE ERRORS FOR DIFFERENT OBSERVERS FOR THE WORST INITIAL CONDITIONS 



B. A REDUCED-ORDER OBSERVER 



The purpose of this example is to illustrate the appli- 
cation of the algorithm to the design of a reduced-order 
observer for a third-order plant. 

Consider the plant characterized by the transfer function 

Y(S) = 1 

S 3 + 4S 2 + 5S + 2 



This plant can be represented by the state and output 
equations 















— 


x 1 (t ) 




0 10 




x x ( t) 




0 


X 2 (t) 


= 


0 0 1 




x 2 (t) 


+ 


0 


x^(t) 




-2 -5 -4 




x^( t) 




1 , 

















y (t) 



u(t) 



( 82 ) 



[1 0 0 ] 



x x (t) 

x 2 (t) 

x 3 (t) 



( 83 ) 



This plant is a third-order system with a single output so 
a second-order observer with arbitrary eigenvalues can be 
constructed. The output matrix C already has the same form 
as the one developed in the reduced-order observer theory. 
Hence, it is not necessary to change the variables, that is, 
for this case the matrices 






B’ = 



ill 



A' 



A 



B 



The matrices A and B can be partitioned as 




Therefore, the reduced-order observer system matrix is 



P 




- -12 




(8H) 



which has the corresponding characteristic equation 



X 2 + ( i ^+g 1 )X + + g 2 + 5 = 0 (85) 

Suppose one selects the starting observer eigenvalues at 
X = -3, -2. This gives the characteristic equation 



(X + 3) (X + 2) = X 2 + 5X + 6 = 0 



( 86 ) 



Equating the coefficients of the same powers of X in ( 85 ) 
and (86) yields 



G 



1.0 

-3.0 



This is the initial guessed value of the observer gain matrix 
(3. If the weighting matrix Q is chosen arbitrarily as Q = I, 
and equation ( 67 ) is solved for V the result is 



V 



36 + 4g 1 + llg 2 + g 2 

11 40 + 42g 1 + 8g 2 + 2g x g 2 + 8g 1 2 



V 



4 - 5e 1 - 



12 40 + 42 gl + 8g 2 + 2 g;L g 2 + 8 gl 2 



6 + 4g. + g 0 + g n 



40 + 42g 1 + 8g 2 + 2g x g 2 + 8g x 2 



(87) 



Due to the symmetry of V only three equations are needed. 
The proposed algorithm gives the solution for the observer 
gain matrix as 

( 88 ) 

if the sum of the eigenvalues of F is limited to -30.0. 

This gain matrix G gives the observer eigenvalues 



G = 



26.0 

-1.73^38 



X 



-25.851 

-4.149 



(89) 



To compare this optimally designed reduced-order 
observer with others having different observer eigenvalues, 
consider two observers with the eigenvalues 



43 



X 



(90) 









-15.0 


, A = 


-20.0 


-15.0 




O 

o 
1 — 1 

I 


— _ 







respectively. The corresponding observer gain matrices are 



«. — 






26.0 




26.0 




^ 0 . ~ 




116.0 




91.0 


- „ 







respectively. In Fig. 9 the actual and estimated states 
using the optimal gains are plotted for the best initial 
error conditions and in Fig. 10 for the worst initial error 
conditions. It is assumed that the initial error condition 
lies on the unit circle. The estimate errors for the three 
observers are plotted on one graph in Figs. 11 and 12 for 
the best and worst initial error conditions respectively. 
Table 3 shows the values of the performance measure for 
observers with different eigenvalues for the worst initial 
conditions. 

C. OBSERVATIONS 

The examples indicate that the observers designed by 
using the proposed technique avoid large errors at a cost of 
having relatively small errors that slowly approach zero. 
This is to be anticipated from the form of the performance 
measure selected. 

The design approach is, however, conservative in the 
sense that the resulting observers function effectively for 
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FIGURE 9 ACTUAL STATES AND ESTIMATES USING THE OPTIMAL GAINS FOR THE BEST INITIAL 
CONDITION 





it 6 



(a) x(2) AND ITS ESTIMATE (b) x(3) AND ITS ESTIMATE 

FIGURE 10. ACTUAL STATES AND ESTIMATES USING THE OPTIMAL GAINS FOR THE WORST INITIAL 
CONDITION 
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(a) ESTIMATE ERROR FOR X(2) (b) ESTIMATE ERROR FOR X (3 ) 

FIGURE 11. ESTIMATE ERRORS FOR DIFFERENT OBSERVER EIGENVALUES FOR THE BEST INITIAL 
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Observer Gains 


Worst Initial 


Values of Perform- 


(Eigenvalues ) 


Conditions 


ance Measure 


&, = 26.0 n = - 25 . 851 \ 




0.122 






1 f 1 


e(0) = 






0.123 


^ « -1.73438 \^ 2 = -H.iHsJ 




- 0.993 






&, = 26.0 A = -15.0 \ 




0.98 






) 


e(0) = 






1.150 


= 116.0 U 2 = -15.0 J 




-0.2 






g, = 26.0 f\ = -20.0 \ 




0.966 








e(0) = 






0.842 


= 91.0 = -10.0 J 




-0.26 

















TABLE 3 . PERFORMANCE VALUES FOR DIFFERENT OBSERVER EIGENVALUES 



situations in which the initial conditions are worst. In 



exchange for this characteristic, performance for less 
difficult initial conditions is compromised. 
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V. CONCLUSIONS 



In this thesis a new method for designing observers has 
been developed. The optimal observer gains are determined 
by solving a minimax problem. The proposed, algorithm is 
easily suited to digital computation and especially appli- 
cable to higher-order systems. 

The observers computed by using the algorithm have a 
tendency to avoid large errors at the expense of small 
errors which slowly decay. As anticipated, the observers 
designed by the proposed technique are efficient for the 
worst initial error conditions . 
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:Cr> = :(r 



fr\ 






:C) = Z 



StET' 
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C V - v V * .■* 



and to the initial state by 



e(t f ) = £(t f ) e ( 0 ) 

Differentiating (A. 4) with respect to t gives 
&(t f - t) = - £(t f - t) F 

for all t^ and t with the boundary condition 
&(tf - t f ) = £(0) = I 



The performance measure for (A.l) is given by 



J = / [ e T (x) Q £(x) ] dx 

0 



Substituting for £(x) using (A. 2) gives 



J = / [ e T (0) <j>. T (x) Q £(x) e(0) ] dr 

0 



Rearranging (A. 8), one obtains 

J = e T (0) [ / /(x) Q £(x) dx ] e_( 0) 

0 

Defining 



V(t) £ / £ T (x - t) Q £(x - t) dx 



then 



V(0) = / £ T (x) Q i(x) dx 

0 



(A. 5) 



(A. 6) 



(A. 7) 



(A. 8) 



(A. 9) 



(A. 10) 



(A. 11) 
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Therefore, we have the final result (29) 



J = e (0) V(0) e ( 0 ) 



( 29 ) 



Differentiating (A. 10) with respect to t by using Leibnitz’s 
rule 



V(t) = 



C £ T (t 



- t) Q £(t “ t) + ± (x - t) Q £(x - t) ] dt 



- ( t - t) Q^(t - t) 



1 • m 

/ C 1 T (t 

t 



, T 



- t) Q £(t - t) + £ (x - t) Q £(x - t) ] dr 



- Q 



(A. 12) 



Substituting (A. 6) for ^(x - t) gives 



V(t) = / ( - [f 1 

t l 



T ^ T ( T - t ) Q £(x - t) ] 



- [± T (x -t)Q£(x-t)F]j dx - Q (A. 13) 



Rearranging (A. 13) gives 



V(t) = - F T [ / 4. T (x - t) Q ±(t - t) dx] 



- [ / £ (x - t) Q £(x - t) dx] F - Q 

t 



(A. 14) 
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Substituting for V(t) using (A. 10) gives equation (30) 



V(t) = - F T V(t) - V(t) F - Q 

Letting t = t„ in (A. 10) gives equation (31) 




Taking the transpose of (A. 10), the right side does not 
change, hence V(t) is symmetric. 



( 30 ) 



(3D 
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